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1 Introduction 

The answers to fundamental science ques- 
tions in astrophysics, ranging from the his- 
tory of the expansion of the universe to the 
sizes of nearby stars, hinge on our ability 
to make precise measurements of diverse as- 
tronomical objects. As our knowledge of the 
underlying physics of objects improves along 
with advances in detectors and instrumenta- 
tion, the limits on our capability to extract 
science from measurements is set, not by our 
lack of understanding of the nature of these 
objects, but rather by the most mundane of 
all issues: the precision with which we can 
calibrate observations in physical units. 

In principle, photometric calibration is a 
solved problem - laboratory reference stan- 
dards such as blackbody furnaces achieve 
precisions well in excess of those needed for 
astrophysics. In practice, however, transfer- 
ring the calibration from these laboratory 
standards to astronomical objects of inter- 
est is far from trivial - the transfer must 
reach outside the atmosphere, extend over 
An steradians of sky, cover a wide range 
of wavelengths, and span an enormous dy- 
namic range in intensity. 

Virtually all spectrophotometric observa- 
tions today are calibrated against one or 
more stellar reference sources, such as Vega, 
which are themselves tied back to laboratory 
standards in a variety of ways. This system's 
accuracy is not uniform. Selected regions of 
the electromagnetic spectrum are calibrated 
extremely well, but discontinuities of a few 
percent still exist, e.g., between the optical 
and infrared. Independently, model stellar 
atmospheres are used to calibrate the spec- 
tra of selected white dwarf stars, e.g. the 
HST system, but the ultimate accuracy of 
this system should be verified against lab- 
oratory sources. Our traditional standard 
star systems, while sufficient until now, need 
to be improved and extended in order to 
serve future astrophsyics experiments. 

This white paper calls for a program to 
improve upon and expand the current net- 
works of spectrophotometrically calibrated 
stars to provide precise calibration with an 
accuracy of equal to and better than 1% 
in the ultraviolet, visible and near-infrared 
portions of the spectrum, with excellent sky 
coverage and large dynamic range. 



2 Science Requiring Precision Cali- 
bration 

The following sections present four science 
investigations that already are or soon will 
be limited by the accuracy of photometric 
calibration. This list is not intended to be 
exhaustive. 

2.1 Expansion history of the Uni- 
verse using Type la supernovae 

In 1998 we learned that the expansion of 
the universe is accelerating, implying the ex- 
istence of a new component of the universe 
dubbed "dark energy". Precise measure- 
ment of the history of expansion and thus 
the properties of dark energy is a major sci- 
ence goal of the next deca de. The Dark 
Ener gy Task Force (DETF) (lAlbrecht et al. 



20061 ) has identified Type la supernovae as 
being one of four principal methods for prob- 
ing the expansion history. 
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Figure 1: Differential magnitude-redshift dia- 
gram for dark energy models with $7, wo, and 
w' = xw a . The difference between models is of 
orde r 0.02 magnitudes (or rou ghly 2%). Models 
from IHuterer &: Linden 120031 . 

Type la supernovae are thought to be 
"standardizable candles" - from observa- 
tions of light curves and spectra, one can 
derive the luminosity of a supernova that is 
the same on average with a scatter of ~15% 
for a single object. Cosmological and dark- 
energy parameters are determined from the 
shape, not the absolute normalization, of the 
Hubble brightness-redshift relationship. For 
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each supernova, its rest-frame B-band flux 
is plotted against its redshift, z. Since the 
rest-frame B-band is seen in different bands 
at different redshifts, the relative zero-points 
of all bands from 0.35 /im to 1.7 /im must be 
cross-calibrated to trace the supernova from 
z = to z = 1.7. 

Planned dedicated experiments, includ- 
ing Pan-STAR RS3, the Dar k Energ y Sur- 
vey (DES) flAbbott et all l200l . the 



sion of a JDEM mission, to a range of possi- 



Large Synop t ic Su rvey Telescope (LSST) 
(llvezic et all 120081 ). and the Joint Dark 
Energy Mission^ (JDEM) and current and 
future observing programs using multipur- 
pose facilities such as the Supernova L egacy 
Survey on CFHT flAstier et all I2006T ). SN 
programs u s ing H ubble Space Telescope 
( iRiess et all 120071 ) a nd James Webb Space 
Telescope (JWST) (IGardner et all 12001 



are or will be focused on collecting accurate 
data for large numbers of supernovae, even- 
tually leading to a data set containing thou- 
sands of objects ranging in redshift from 
to 1.7. 
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Figure 2: Simulated SNIa data from one ver- 
sion of a JDEM mission compared with pre- 
dictions from a range of Dark Energy Models 
(Derived from lWeller fc Albrechtl f|200H )). 

The power of using SNela out to z~ 1.7 
for measuring the cosmological parameters 
is demonstrated in Figure [21 which compares 
the expected (simulated) results of one ver- 



1 http://pan-starrs.ifa.hawaii.edu/public/ 
2 http:/ /jdem. gsfc.nasa.gov/ 



ble d ark energy models flWeller fc Albrecht 



20011 ) . This calculation is based on 2000 



SNe la measured in the range 0.1 < z < 1.7 
plus 300 low-redshift SNe l a from, e.g., the 
Near by Supernova Factory flAldering et al , 
120021 ). The simulated data have a statisti- 
cal accuracy that is capable of distinguishing 
models whose predictions differ by as little 
as 2% over the full range of redshifts. 

However, to make full use of the data, 
systematic errors must be comparable to 
or smaller than the statistical errors. The 
NASA-DOE Joint Dark Energy Mission's 
Reference Mission specifies that, over the 
fullwavelength range of 0.35 < A < 1.7 /im, 
a photometric uncertainty of 0.5% per oc- 
tave is required for the mission to reach its 
target Figure of Merit. Achieving this level 
of precision at the (faint) flux levels of the 
redshifted SNe requires a transfer of the ab- 
solute calibration from bright standard stars 
to fainter calibration standard stars which 
can be directly observed by the DE missions. 

2.2 Growth Of Structure 
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Figure 3: A figure demonstrating the technique 
of photometric redshifts. For each galaxy type, 
the dashed line joins together points that mark 
the colors at a particular redshift. From lower 
right to upper left, the redshift increases from 
to 0.6 (Eisenstein et al. 2001). 

A potentially powerful technique for mea- 
suring the growth of structure in the uni- 
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verse is to use gravitational weak lens- 
ing combined with photometric redshifts of 
galaxies to study the statistical properties of 
the mass distribution clS cL function of red- 
shift. The history of growth of structure 
provides another approach to measuring the 
properties of Dark Energy. Current and 
future experiments that propose to collect 
data for such studies include Pan-STARRS, 
DES, LSST, and JDEM. 

A simplified description of this approach 
is as follows. One identifies a set of galax- 
ies at the same approximate redshift and 
measures the distortions in the shapes of 
these galaxies induced by gravitational lens- 
ing from the intervening mass distribution 
(such as clustering). A single set of galax- 
ies measures the properties of the integrated 
mass distribution along a line of sight. By 
selecting a second set of galaxies at, e.g., a 
higher redshift, and measuring the changes 
in lensing-induced shapes relative to the first 
set, one obtains information about mass 
structures in a slice of space between the two 
sets of galaxies. Thus, in a process analo- 
gous to tomography, one can build up a view 
of the mass structures and how they change 
clS 8b function of redshift. 

A key necessity in this approach is the 
use of multicolor, broadband photometry of 
galaxies as a "low-resolution spectrograph" 
to estimate redshifts (Fig. [S]). Because 
the intrinsic spectral energy distribution of 
any galaxy is not known a priori, one must 
rely on matching a set of redshifted tem- 
plate spectra to the measured photometry 
of a galaxy and utilizing a "training set" of 
galaxies with known redshifts to calibrate 
the templates. 

Spectrophotometric calibrations are used 
to convert the template spectra to pre- 
dictions of galaxy magnitudes and colors. 
Ideally, the training set would span all of 
parameter space, but in reality there will 
always be galaxies that can be measured 
photometrically but are too faint to mea- 
sure spectroscopically. Accordingly, the 
LSST project has developed a two-pronged 
approach to obtain photo metric redshifts 
from its multicolor data set (jConnolly et al . 
120061 ). and established a requirement on 
spectrophotometric calibration of 1% (1.5% 
in the UV), with design goals that are twice 



as goocQ. 

2.3 Stellar Populations In Elliptical 

Galaxies 



Although elliptical and SO galaxies are 
only a small fraction of all galaxies, they are 
notable for having very similar stellar pop- 
ulations, as reflected in their uniformity of 
colors. With the advent of large, multicolor 
surveys using digital detectors, these objects 
can be identified over a range in redshift 
and used for cosmological studies. Thus, the 
red galaxy spectroscopic sample in SDSS has 
been used to detect acoustic baryon oscilla- 
tions (lEisenstein et a/I 120051 ). Additionally, 
optical detection and measurement of galaxy 
clusters has seen a resurgence of interest 
due to the ability to identify galaxy clusters 
based on the "red sequence" of these types of 
galaxies. In low redshift clusters, the colors 
of early-type galaxies are remarkably uni- 
form, showing a scatter of ju st 5% in color s 
such as SDSS g — r and r — i dKoester et al 



255% . The SpARCS survey (IWUson et al 



20081 ) has shown that clusters with similar 
galaxy content exist out to redshifts of at 
least 1.34. Galaxy clusters will be detected 
and measured by nearly every current and 
future imaging survey conducted for weak 
lensing. Galaxy cluster counts have been 
identified as a third method for measuring 
dark energy by the DETF. 




Figure 4: Black squares: colors of ellipti- 
cal galaxies as a function of redshift; Red 
crosses: Passively evolving stellar population 
model (Eisenstein et al. 2001). 



3 www. lsst.org/Science/docs/SRD.pdf 
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Since a cluster has anywhere from 10 to 
100 members, the mean color of galaxies can 
be measured with extremely high precision. 
By comparing the colors over a range in red- 
shifts, it should be possible to make accu- 
rate models of the stellar populations and 
infer their evolution over a significant frac- 
tion of the age of the universe. The limit 
on the accuracy of these models will be set 
by the ability to self-consistently calibrate 
the galaxy photometry over the optical and 
near-IR bands. Conceivably one could take 
advantage of data calibrated at better than 
1% accuracy. Figure |4] demonstrates the 
precision with which elliptical galaxy colors 
can be measured and compared with stellar 
synthesis models. 
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Figure 5: Color-color diagram for DA white 
dwarfs. Open and filled circles are observations 
of stars with measured distances. Solid lines 
are predictions from a grid of models with con- 
stant gravity or constan t effec tive temperature 
(IHolberg. J fe Bergeron! . HH). 



2.4 Stellar Structure 

The fundamental parameters of stars, in- 
cluding mass, radius, metallicity, and age, 
are inferred by matching accurate models 



of stellar atmospheres to calibrated spec- 
troscopic data and thus determining effec- 
tive temperature, surface gravity, composi- 
tion and, if necessary, interstellar redden- 
ing. For stars with relatively simple atmo- 
spheres such as hydrogen white dwarfs, at- 
mosphere models are thought to be quite ac- 
curate and can be used to predict photomet- 
ric parameters (Fig. [5]) and, in combination 
with stellar interior models, the radii and 
absolute luminosities as well. By combining 
these data with photometric measurements, 
it is possible to predict distances. A com- 
parison of these predictions with measured 
trignometric parallaxes for those stars with 
such measurements s hows excellent agree- 
ment (IHolberg et all 120081 ). If calibrations 



can be improved to the level of 1% and with 
more stars (such as will be measured with 
GAIA), it will be possible to make mean- 
ingful tests of 3-D spherical models, derive 
masses directly, and make more quantitative 
tests of evolutionary models. 

3 Flux Calibration & Standardiza- 
tion 

Ultimately, observed astrophysical fluxes 
must be converted to physical units. Three 
of the most common methods of determining 
the absolute fluxes are through comparison 
to standard stars (e.g. solar analog stars), 
stellar atmosphere models, and certified lab- 
oratory standards. But, the existing preci- 
sion of each of these methods is inadequate 
to meet the requirements of the science de- 
scribed in the previous section. 

3.1 Solar Analog Stars 

Use of solar analog stars standard 
source relies upon the star having the same 
intrinsic SED as the sun. Unfortunately, no 
star is a true solar analog. Even G-type stars 
with the most-closely matching visible spec- 
tra can differ by a few percent. In addition, 
uncertainties in the solar SED itself are 2- 
2003h . 



3% flThuillier et al 



3.2 Stellar Atmosphere Models 

UV and visible astrophysical fluxes are of- 
ten normalized to an absolute flux using a 
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set of hot, white dwarfs (WDs) whose mod- 
els are tied to Vega's absolute flux at 5500 
Angstroms, as determined through direct 
comparison to a black body reference. 

Stellar atmosphere models are currently 
the preferred method for calibrating stel- 
lar fluxes due to the agreement between 
the models and the observations as well as 
the increased resolution of both the models 
and the data. Use of these pure hydrogen 
WD stars simplifies the computation and 
improves the precision by eliminating one of 
the most difficult steps in atmospheric mod- 
eling - that of including the blanketing from 
the plethora of metal lines. 

To obtain the absolute flux and its un- 
certainty for an unreddened WD, medium- 
resolution high S/N (>50) observations of 
the Balmer lines are fit to model hydrogen 
line profiles to determine the effective tem- 
perature, the gr avity, and the associ ated un- 



certainties (e.g.. iFinley et all 119971 ). Then 



the best-fit model and the models at the ex- 
tremes of the uncertainty in T e fj and logg 
determine the nominal flux and uncertainty 
in the shape of the flux distribution. These 
model fluxes are normalized to V-band Lan- 
dolt photometry. 

The three primary WD standards of HST 
CALSPEC network are internally consis- 
tent to an uncertainty level of 0.5% in the 
visible with localized deviations from mod- 
els rising to ~1% over the 4200-4700 A 
spectral range, and a ± 1% uncertainty i n 

El ' 



Bohlinl 120071) 



the NIR (1-2 //m) (Fig. 
Current uncertainties in the extensive NIR 
(1.0 < A < 1.7 fjm) network of standa rd 



stars are 



e.g. 



Cohen et al\ Il992al lbl. 



20031 : ICohenll2007l h 

Any systematic modeling errors that 
equally affect the shape of the flux distribu- 
tions of all three WD stars cannot be ruled 
out and would make the actual error larger. 
Differences between the continua of the LTE 
and NLTE models place a lower limit of 2% 
on the uncertainty in the 0.35—1.7 /im range 



for these standards. 

In the NIR, astrophysical fluxes are often 
normalized to A-star models, where the ac- 
curacy of the best A-star models rivals that 
of the pure hydrogen WD models. Abso- 
lute photometry of Vega is used to normalize 
the SEDs of these stars to an absolute flux 



scale. iRieke et al\ (120081 ) tested the agree- 
ment of IR standard star calibrations and 
models based on direct absolute measure- 
ments of A0V stars versus the sun and exam- 
ined the impact of extrapolating the IR data 
into the visible. The data were found to be 
consistent, permitting flux calibrations with 
an accuracy of ~2% between 1 and 25 /im. 



VEGA: STIS-black, Kuruez(9550K)-red, Kurucz(9400K)-green 
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Figure 6: Uncertainties in the absolute flux 
for Vega: HST/STIS observations (black line: 
IBohliriEnnftlBohlin fc (hffiBl2TO . the Ku- 
rucz stellar model with Tg/-j=9400K (green), 
and the Kurucz stellar model at 9550 K (red) 
are compared. The observations exhibit better 
agreement with the cooler model at the longer 
and shorter wavelengths. The hotter model 
agrees better with the measured flux by ~ 1% 
at 4200-4700 A. 

3.3 Certified Laboratory Standards 

Photometry of Vega has been absolutely 
calibrated against terrestrial observations of 
certified laboratory standards (e.g. tungsten 
strip lamps, melting point black bodies) to 
provide the normalization for the network of 
stellar models and templates that are used 
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as practical absolute standards. These ab- 
solute calibrations to standard sources were 
difficult and subject to large systematic un- 
certainties due primarily to the large and 
variable atmospheric opacity. 

Discrepancies of > 10% in Vega's flux 
exist at 0.9— 1/im, whereas the measure- 
ments from 0.5— 0.8 ixm agree to ~ 1% 
f lBohlin fc Gillilandll2004 lHayedll985h . Be- 
yond 1 /im, windows of low water va- 
por absorption have b een used fo r ab- 



Selbv et all 11983 



solute photometry ( e.g 
Mountain et a/.lll985l ). 

Currently, the uncertainty in the standard 
star flux calibration network relative to the 
fundamental laboratory standards exceeds 
1%. 

Certified Detectors: The calibration 
precision of photodetectors has greatly im- 
prov ed since early pioneer i ng measurements 



Oke and Schild 1970l ; Hayes fc Latham 



19751 ). Current NIST ~2cr uncertainties in 
the absolute responsivity of standard detec- 
tors are ~ 0.2% for Si photodiodes and 0.5% 
for N IR photodiodes fjLarason and Houston! . 
20081 ). This increased precision in the pho- 
todetector calibration, ease of use, and re- 
peatability, now make standard detectors 
the calibrator of choice. 

3.4 Extension to Standard Star Net- 
works 

The basic techniques and methodologies 
for extending one fundamental standard 
candle to a network of stellar standards 
are well established. This extensive net- 
work of stellar standards is fundamentally 
tied to the sun or to Vega, e.g SDSS 
successfully established a network of stan- 
dard stars s panning the visible range from 
0.3-1.0 u,m flSmith et all l2002t ) with abso- 
lute fluxes based on Vega using BD+17°4708 
as an interme diate (V=9.5 mag) trans- 



fer standard (IFukugita et all Il996l ). Even 
the ICohen et al\ (|1992al ) absolute standard 
models of Sirius are tied to Vee; the un- 
derlying standard. 



Vega is far too bright to be observed di- 
rectly by the current class of 4-m telescopes 
and even with most 2-m telescopes using 
state-of-the-art instruments. Its use as a 
standard is further complicated by its proto- 
planetary disk which contributes to IR flux 
measurements. In addition, clS db pole-on 
rotator its surface temperature and gravity 
vary dramatically from the p ole to equator 
(e.g. lAufdenberg et all ( 120061 )). This intro- 
duces complexity into accurately and pre- 
cisely representing its flux with robust stel- 
lar atmosphere models. Furthermore, un- 
certainties in atmospheric corrections have 
resulted in wavelength dependent uncertain- 
ties in Vega's intrinsic flux. Thus, Vega is 
not suitable as a modern astrophysical flux 
standard. 

NIST standards have been transferred to 
observations of other stars, but the level of 
uncertainty in the flux measurements have 
precluded their widespread use (e.g. HZ43 
and G191B2B: ~4% precision, iKruk et al 



( 119971 )). An exception to this was the Mid- 
course Space Experiment (MSX) which ob- 
served eight standard stars, including Vega, 
in the infrared and directly tied these ob- 
servations to inflight me asurements of em - 
misive reference spheres (IPrice et all 120041 ) . 
These measurements resulted in corrections 
(Sirius: 1%) and caveats (Vega: flux ex- 
cess). These MSX observations were limited 
to bright, typically Kill and Mill, stars in 
six selected NIR/IR bandpasses. Thus, the 
need for a sample of absolutely calibrated as- 
trophysical standards spanning a broad dy- 
namic range in flux and wavelength (UV 
through NIR) persists. 

Current astrophysical problems need a 
precise (better than 1%) network of astro- 
physical flux standards spanning a wide dy- 
namic range. This enables scientists to take 
advantage of the capabilities of current and 
future telescopes and the instruments that 
were developed to address pressing scien- 
tific questions. New, direct measurements 
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of standard stars tied directly to fundamen- 
tal NIST standards are required. 

3.5 Current Status & Future 
Prospects 

Although the relative photometry of ob- 
jects in a single CCD exposure can be bet- 
ter than 1%, this level of precision is not 
achieved for the relative fluxes of sources 
i n diff erent fields of view. IStubbs fc Tonry 
(120061 ) reviewed systematic uncertainties 
that plague ground-based observations, dis- 
cussed the challenges associated with char- 
acterization of atmospheric transmission 
and the removal of instrument artifacts, 
and presented a method for achieving pho- 
tometry with fractional uncertainties. Us- 
ing precisely calibrated photodiode detec- 
tors in concert with a wav elength tun- 
able l aser illumination source, IStubbs et al. 
(120071 ) demonstrated the success of their 
methodology in measuring the instrument 
transmission and established the capabil- 
ity of standard detectors fundamental 
metrology to achieve precise and accurate 
photometry. 

Other programs are also making con- 
certed efforts to character ize instrument 
perfor mance (e.g ASTRA lAdelman et al. 
(120071 )). however, the need to monitor and 
correct for atmospheric transmission on 
short timescales persists. One approach 
(e.g. Pan-STARRS, LSST) uses a dedi- 
cated telescope to monitor the atmosphere 
throughout the night to enable corrections 
for science observations at the neighboring 
facility. 

Direct, absolute calibrations of stellar 
fluxes measured above the Earth's atmo- 
sphere are also being pursued. A recently 
approved sub-orbital program, ACCESS: 
Absolute Color Calibration Experiment fo r 



20081 . 12007m 



Standard Stars (IKaiser et al. 
will transfer NIST absolute detector stan- 
dards to additional standard stars with bet- 
ter than 1% precision over the -3500A - 
1.7^m bandpass at a spectral resoloving 



power of ~500. However, due to the limited 
time above atmosphere for rocket flights, 
these measurements will be limited to a few 
stars brighter than ~\0 th magnitude. 

The scientific impact of a standard star 
network based on the absolute calibration 
of stars too bright to be observed with the 
premier telescopes needs to be addressed. 
A modern calibration program should ex- 
tend direct flux measurements to fainter 
sources, encompass a broad spectral range 
(UV through the IR), ensure robust results 
through the support of independant calibra- 
tion programs, and provide technology sup- 
port to execute these programs. 

In conclusion, we stress the need for a cal- 
ibration program that supports the science 
of the 21st century. 
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